Mass and redshift dependence of dark halo structure 



D.H. Y.P. Jing^'^, H.J. Mo^ G. BSrner^ 



ABSTRACT 



Using a combination of N-body simulations with different resolutions, we 
study in detail how the concentrations of cold dark matter (CDM) halos depend 
on halo mass at different redshifts. We confirm that halo concentrations at the 
present time depend strongly on halo mass, but our results also show marked 
differences from the predictions of some early empirical models. Our main result 
is that the mass dependence of the concentrations becomes weaker at higher 
redshifts, and at ^ > 3 halos of mass greater than 10^^h~^MQ all have a similar 
median concentration, c ~ 3.5. While the median concentrations of low-mass 
halos grow significantly with time, those of massive halos change only little with 
redshifts. These results are quantitatively in good agreement with the empirical 
model proposed by Zhao et al. which shows that halos in the early fast accretion 
phase all have similar concentrations. 

Subject headings: cosmology: miscellaneous — galaxies: clusters: general — 
methods: numerical 



High-resolution A^-body simulations have shown that the density profiles of cold dark 
matter (CDM) halos can be described reasonably well by a universal form. 



where is a characteristic "inner" radius, and ps is the density at (Navarro, Frenk and 
White 1996, 1997; NFW hereafter), although there is still debate about the exact value of 
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the inner slope (e.g. Pukushige &: Makino 1997, 2003; Moore et al. 1998; Jing & Suto 
2000). For a halo of radius R, this profile is characterized by the concentration parameter 
c = R/vg, which is found to be dependent on halo mass (smaller halos have, on average, higher 
concentrations). NFW developed a simple model to account for this mass dependence. The 
NEW model for the mass dependence of c was later found to be inconsistent with results 
obtained from a simulation of the concordance low-density CDM model (LCDM) at high 
redshift (Bullock et al., BOl hereafter). BOl found that c is proportional to the cosmic scale 
factor a for a given halo mass. Based on this, BOl and Eke et al. (2001; EOl hereafter) 
proposed new empirical prescriptions to predict c as a function of redshift and halo mass. 

It must be noticed, however, that all these prescriptions only give the mean concentra- 
tion of all halos of a given mass at a given redshift. Since the density profiles are found to 
vary significantly from one halo to another even for a given halo mass (Jing 2000; BOl), it is 
important to have a recipe to predict c for individual halos. Jing has examined the density 
profiles for halos in different dynamical states, and found that the halo density profiles are 
closely related to the halo formation history. The connection between halo concentration 
and halo formation history was explored further by Wechsler ct al. (2002; W02). Assuming 
that c = 4.1 at their defined "formation redshift" and c oc a, W02 found that their model 
prediction for c is in good agreement with their simulation results for the LCDM model. 

In a recent paper, Zhao et al. (2003; hereafter ZMJB) found that, for a given halo, 
there is a tight correlation between the inner scale radius and the mass within it, Mg, for 
all its main progenitors, and that this correlation can be used to predict the concentration 
of a dark halo at any time without making any ad hoc assumption about the form of the 
mass accretion history. The ZMJB model predicts that the evolution of c of individual halos 
are not just a function of a (such as c oc a as W02 assumed), but tightly connected to their 
mass growth rate: the faster the mass grow, the slower the c increase. 

In this Letter we use a combination of high-resolution simulations of different boxsizes to 
directly explore the halo structures for a wide range of halo mass in a wide range of redshifts. 
This allows us to study in detail how halo concentration depends on halo mass at various 
redshifts, and to test the accuracy of the various empirical models mentioned above. We will 
show that at high redshift the simulated mass dependence of c is much different from some 
previous results, and is quantitatively in good agreement with ZMJB prediction. 
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2. Simulation results 

The cosmological simulations used in this paper are generated with a parallel- vectorized 
Particle-Particle/Particle-Mesh code (see Jing & Suto 2002). The concordance CDM model 
with the density parameter Qq = 0.3 and the cosmological constant Aq = 0.7 is considered. 
The linear power spectrum has the shape parameter F = D,oh — 0.20 and the amplitude 
as — 0.9, where h is the Hubble constant in 100kms~^Mpc~^, and as is the rms top-hat 
density fluctuation within a sphere of radius 8/i~^Mpc at the present time. We use 256^ 
particles for the simulation of boxsize L = 25/i~^Mpc, and 512^ particles for the other two 
simulations of L = 100 and 300^~^Mpc (Table 1). The simulations with L = 25 /i~^Mpc 
and 100 h^^Mpc have been evolved by 5000 time steps with a force softening length r] (the 
diameter of the S2 shaped particles, Hockney & Eastwood 1981) equal to 2.5/i^^kpc and 
10/i~^kpc, respectively. As a result, there are many halos with more than 3000 particles 
in these simulations, and these halos are resolved similarly to or better than the individual 
halo simulations in early studies of the density proflles (e.g. NFW). It has been shown that 
the resolution at this level is sufficient for determining the concentration parameter (Jing 
2000; EOl), though it may not be good enough for addressing the issues with regard to 
the slope of the density profile in the central region of a halo (e.g. Fukushige & Makino 
1997, 2003; Moore et al. 1998; Jing & Suto 2000). The simulation of L = 300 /i^^Mpc is a 
typical cosmological simulation, evolved by 1200 steps and with a force softening length of 
30/i^^kpc. The halo sample constructed from these simulations is big, which is essential for 
accurately determining the mean halo concentration. There is a sufficient overlap in mass 
between the halos of different simulations, from which the resolution effect can be rehably 
estimated. The halos are defined according to the spherical virialization criterion (Kitayama 
& Suto 1996; Bryan & Norman 1998), so the radius of a halo in this paper is the virial 
radius rvir- The halos are identified from simulations using the potential minimum method 
as described in Jing & Suto (2002), and the particle with the minimum potential in each 
halo is chosen as the halo center. We use all halos identified this way without applying any 
further selection criteria. 

Table 1: A summary of simulation parameters 



simulation 




box size 




rj 


-^initial 






/i-^Mpc 


h-^MQ 


h~^kpc 




LCDM025 


256=^ 


25 


8.0 X 10^ 


2.5 


72 


LCDMIOO 


5123 


100 


6.4 X 10^ 


10 


72 


LCDM300 


512'^ 


300 


1.7 X IQi" 


30 


36 



-4- 



2.1. The mass accretion history 

Following ZMJB, we construct the main branch of the merger tree for each halo identified 
at redshift z = 0, and work out the mass accretion history. In Figure 1, we plot the 
mass accretion histories for 20 randomly selected halos at each of the following mass scales: 
Mvir,o = 7 X lOio/i-^M©, 1.3 X lO^^/i-iMo, and 1.4 x IO^^H-^Mq. As in ZMJB, the mass 
accretion history of each halo is divided into a fast accretion phase and a slow accretion phase, 
and we denote the transition redshift between the two by Ztp. With the mass Mvir(-2) in units 
of Mvir,tp = M^hiztp) and the physical virial density pviv{z) in units of 
we found that the mass accretion history has a universal form for halos of different mass, 
and the mean accretion history can be accurately represented by the thick smooth line for 
the halo masses covered by the simulations (from 7 x IO^^H'^Mq to 1.3 x IO^^H'-^Mq). 
Mathematically, the average accretion history can be expressed in the form: 

where x = Pvir,tp/Pvir(^) and a = 0.75 (0.42) for the fast (slow) accretion phase. The universal 
mass history has interesting implications for galaxy formation and halo structure formation. 
Here because of the limited space, we will only discuss its implications for halo structure in 
^3. 



2.2. Halo concentrations 

We select halos with more than 800 particles at redshifts 2 = 0, 1, 2 and 4, and 
determine the concentration parameter for each of them by fitting the density profile to 
the NFW form (see ZMJB for the fitting procedure). The halos are grouped in mass bins 
of Alogio^vir = 0.2, and the median concentration is calculated for each mass bin. The 
median concentrations determined this way are presented in Figure 2 for masses larger than 
6.4 X lO^°/i~^M0, together with their errorbars (standard deviation among different halos 
in the bin divided by the square root of the halo number). Note that there is always quite 
a large overlap in halo mass between simulations of different boxsizes, and that the median 
concentrations from different simulations are in good agreement in the overlapping mass 



^For the LCDM model, /Ovir(-z) is 180 times at high redshift and 101 times at z = of the critical density. 
Both pvir{z) and 1 + z can be used to denote the cosmic time for a given cosmology. We found the relation 
between Mvir and pvir{z) is better behaved than the Afvir - z relation, especially for low-density universes at 
z ~ 0. 



- 5 - 



F 1 I I I I I 



• 1—1 

> 



10 E 



1 - 



0.1 E 



0.01 E 



0.001 



— I — rrm 

M. 



I I 1 1 1 1 



I I 1 1 1 1 



"1 — rrm 



I M 1 1 im 



M^.^o~7.0xlOio h-^M. 




// a = 0.75 (x<l); a = 0.42 (x>l) 
predicted 



I i/i 1 1 1 























0.001 0.01 



0.1 

X = p 



10 



100 1000 



vir.tp 



p 



vir 



Fig. 1. — The mass accretion history of dark matter halos along the main branch. Twenty 
halos are randomly selected from the simulations at each mass indicated at the top. We use 
the physical density within the virialized halo pv\r{z) as the time variable. Both the halo 
mass Mvir(-2) and Pm{z) are scaled to their quantities at the turning point, Mvir.tp and Pvir,tp- 
The scaled mass accretion history on average is well represented by Eq.(2) independent of 
the halo mass. The thick smooth dashed line is the prediction of ZMJB for the concentration 
Cvir according to the mean accretion history (the right vertical coordinate for Cvir)- 
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range. This agreement suggests that our results are not significantly affected by the finite 
numerical resolution, because for a given halo mass both the mass and force resolutions 
decrease with the increase of the simulation boxsize. 

As one can see, at low redshift, the median concentration decreases rapidly with halo 
mass, from c ~ 20 for Mvir ~ lO^^h'^MQ to c ~ 4 for Mvir ~ lO^^h'^MQ. This mass 
dependence is similar to that found in earlier analyses (e.g. NFW; Jing &: Suto 2000; BOl 
). The mass dependence is weaker for halos at higher redshifts and becomes insignificant at 
z > 3. This change in behavior is mainly due to the fact that the median concentrations 
of small halos decrease rapidly with increasing redshift while the concentrations of massive 
halos change little. Note that the decrease of c with z is slower at higher z and there seems 
to be a minimum value c ~ 3 for the median concentration of dark halos. This is true even 
for a few halos at 2; ~ 9 that are not included here. 



3. Compcirison with empirical models 

Based on results from numerical simulations, ZMJB found that the scale radius of a 
halo and its scale mass Ms (i.e. the mass within r^) are tightly correlated, with a relation 
well represented by a simple power law: 

where Mg^ and rs,o are the scale mass and scale radius at some chosen epoch. The value of 
a is found to be 0.52 in the slow accretion phase, and 0.64 for the rapid accretion phase. As 
shown in ZMJB, this Mg-Vg relation can be used to derive c from the halo mass accretion 
history according to 

[lu(l + r) -rV(l + r-)]r--^° 
[ln(l + co) -Co/(l + Co)]co^ 

This relation can be calibrated by fixing ctp at 2; = Ztp. ZMBJ have calibrated Ctp with 
five high resolution halos, and adopted Ctp — 4.0. With our current large sample, we find 
Ctp = 3.5 to be more accurate. 

With the universal mass accretion history obtained from our simulations, this recipe 
can be used to predict the median concentration as a function of redshift. The result is 
shown as the smooth dashed line in Figure 1. As one can see, the halo concentration has 
a value about 3.5 in the fast accretion phase, and scales roughly as c oc p~^^^(z) in the 
slow-accretion phase. This is consistent with our simulation results, that halos with masses 
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Fig. 2. — The median concentration of halos as a function of the halo mass. Errorbars are 
the standard deviation among halos of the same mass devided by square root of the halo 
number. The black symbols are for the results measured in the simulations, the blue and red 
lines are the predictions of ZMJB using the mass accretion histories from the simulations 
(blue) and from the PINOCCHIO model (red), and the dashed and dotted lines are the 
predictions of the models of Bullock et al. and Eke et al., respectively. For comparison, the 
simulation results of Bullock et al. for z = 4 are plotted in the lower right panel as green 
pentagons. 
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between 10^^ - lO^^/i~^M0 have median concentrations independent of 2; at 2; > 1. Most 
of these massive halos are in the fast accretion phase at these high rcdshifts. Note that 
the median concentration obtained from the simulations never drops below 3 even for the 
most massive halos at the highest redshift probed by our simulations, in agreement with the 
ZMJB model. The strong increase of c with decreasing z for low-mass halos at low redshift 
observed in the simulations is also consistent with the model prediction, because most of 
those halos are in their slow accretion phases. 

We generate samples of mass accretion histories from the simulations, and apply the 
ZMJB model to predict the concentrations for each of these halos (Figure 2). The model 
prediction reproduces well the mass and redshift dependence of halo concentrations obtained 
from the simulations. The distribution of the concentrations for given halo mass and redshift 
is well described by the log-normal distribution with crinc ~ 0.3 (Jing 2000, BOl). We also 
compute mass accretion histories using the PINOCCHIO code of Monaco et al. (2002). 
This code identifies dark matter halos and their merging histories by applying an ellipsoidal 
collapse model to an initial cosmic density field. It has been shown to be quite accurate in 
reproducing many properties of the halo population. All parameters are kept the same as 
in the simulaitons, and also we apply the ZMJB model to predict the halo concentrations. 
Again the agreement with our simulation results is satisfactory with an accuracy better than 
10% (Figure 2). In a forthcoming paper (Zhao et al., in preparation), we will show that 
the ZMJB prediction is valid for a wide range of cosmological models, including the SCDM 
model, an OCDM model, and scale-free models. 

The increase of the halo concentration with decreasing redshift in the slow accretion 
phase is qualitatively consistent with the relation c{My\^^z) oc (1 + z)~^ found by BOl, 
because p^,\^{z) is approximately proportional to (1 -|- zf" . It is however important to note 
that the evolution we obtained is along the main branches of merger trees, while the relation 
c(Mvii., z) oc (1 -|- z)"-*^ obtained by BOl is for the median concentration of a halo population 
of a given mass Mvu- There is a marked difference of the prediction of the ZMJB model 
from that of BOl: while ZMJB predicts that c does not change in the fast accretion, BOl 
predicts c oc (1 -|- z)~^ for a given mass. 

In Figure 2, we compare our results to the predictions of the empirical models given by 
BOl and EOl. The BOl model agrees with our simulation results well at redshift z — ^ for 
Mvir < lO^'^/i^^M©. Note that this is approximately the mass range that their simulation 
can effectively explore; their simulation uses 256^ particles in a box of 60/i^^Mpc. The BOl 
model also agrees with the redshift dependence of c for low-mass halos, but it fails to match 
our simulation results at the high mass ends. The model of EOl fits our simulation data 
better for 2; = 1 and 2, but worse for 2; = than the BOl model. Since the Eke et al. model 
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adopted a redshift-dependence of c similar to that of BOl model, it also underestimates the 
concentration for massive halos at high redshift. 

It should be pointed out that the models of BOl and EOl both match their own sim- 
ulations to redshift 4 and 2, respectively, and so there seems to be a discrepancy among 
the different simulation results. Our results are consistent with the simlation results of EOl 
over the mass and redshift ranges probed by their simulations. Note that the very low 
concentrations obtained by EOl are for a warm dark matter spectrum. There is a marked 
difference between BOl's simulation results and ours for z > 2; as comparison we plot in 
the low-right panel of Figure 2 their simulation results for 2; = 4 (green pentagons). As 
one can see, the discrepancy between BOl and our simulation results becomes significant 
at Myir ^ 5 X 10^^h~^MQ for this redshift. Unfortunately, the origin of this discrepancy is 
unknown. Our 100/i~^Mpc box simulation has a mass resolution slightly better than their 
main (60/i~^Mpc) simulation, and in terms of halo number, our sample is more than 4 times 
larger. Since the number density of massive halos is low at high z, and since halo concen- 
tration can differ substantially for halos with the same mass at the same redshift, a large 
sample might be crucial to get reliable results. We are confident about our results, because 
our halo sample is large and our simulations with different resolutions and boxsizes agree 
with each other very well. 

As mentioned earlier, the consistency of the results obtained for different simulation 
boxsizes indicates that our determination of halo concentrations should not be affected sig- 
nificantly by the limited simulation resolution. As shown by Moore et al. (1998) and Die- 
mand et al. (2003), the finite resolution should reduce the halo concentration. Comparing 
our results for different boxsizes, it appears that the resolution effect can lead to an under- 
estimate of c by < 5% at the lower halo mass end in each simulation. The slight systematic 
difference between the simulation results and the ZMJB model predictions may be caused 
by this resolution effect, and correcting for this may give a better agreement between the 
simulation data and model predictions. We have also examined the validity of the NEW 
profile for halos that are in the fast accretion phase, and found that most of these halos can 
be fitted by this profile (Zhao et al. in preparation). 

4. Discussion and conclusions 

We have studied the dependence on mass and redshift of the concentration of cold dark 
matter (CDM) halos in high resohition simulations, and discovered that at early times the 
mass dependence of halo concentrations which is pronounced at present, becomes insignifi- 
cant, and at z > 3 halos of mass > 10^^h~^MQ have in the mean the same density profile 
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with c ~ 3.5. Our results indicate that the median concentration of halos cannot decrease 
with redshift or /and halo mass to a value less than ~ 3. Massive cluster halos at present 
have higher concentrations than some previous models predicted. 

The good agreement between the results of c from different simulations demonstrates 
that the simulation resolution effect has been well controlled in our analysis. 

For the mass accretion histories of halos with masses from 10^° to IO^^H^^Mq, we have 
found that they are well expressed by the universal function [Eq([2])]. 

All these results can be quantitatively matched by the empirical model of ZMJB. In 
that model the concentration of a halo is related to its mass accretion history through the 
scaling relation they found for and Mg. It predicts that all halos in the fast accretion 
phase have similar concentrations, regardless of their mass or their redshift . The ZMJB 
model reproduces our results both in the fast and the slow accretion phases. Since directly 
modelling halo density profiles in numerical simulations is both expensive and time consum- 
ing, the ZMJB model provides a practically useful technique for modeling internal structures 
of individual CDM halos. 

While the model of Bullock et al. (BOl) agrees with our results for halos in slow 
accretion, it seriously underestimates the concentration of halos in the fast accretion phase. 
The models of EOl have the same weakness as BOl, since they all adopt a similar assumption 
that c(Mvii., a) oc a. 

The implications for galaxy formation models, and for the interpretation of observations, 
such as strong and weak lensing surveys, are interesting. According to our results, massive 
cluster halos at the present time, galactic halos at 2; ~ 3, and halos of the first collapsed 
objects in the universe at z ~ 15 (such as POP III stars) should all have about the same 
concentration. 
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presented in this paper were carried out at ADAC (the Astronomical Data Analysis Cen- 
ter) of the National Astronomical Observatory, Japan. The work is supported in part by 
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